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Following the discovery of Bose-Einstein con- 
densation (BEC) in ultra cold atoms [E. Gosta, 
Nobel Lectures in Physics (2001-2005), World 
Scientific (2008)], there has been a huge exper- 
imental and theoretical push to try and illumi- 
nate a superfluid state of Wannier-Mott excitons. 
Excitons in quantum wells, generated by a laser 
pulse, typically diffuse only a few micrometers 
from the spot they are created. However, Bu- 
tov et al. [1] and Snoke et al. [2] reported lumi- 
nescence from indirect and direct excitons hun- 
dreds of micrometers away from the laser exci- 
tation spot in double and single quantum well 
(QW) structures at low temperatures. This lu- 
minescence appears as a ring around the laser 
spot with the dark region between the spot and 
the ring. Developing the theory of a free super- 
flow of Bose-liquids we show that the macroscopic 
luminesce rings and the dark state observed in 
[1, 2] are signatures of the coherent superflow of 
condensed excitons at temperatures below their 
Berezinskii-Kosterlitz-Thouless (BKT) transition 
temperature [3]. To further verify the dark exci- 
tonic superflow we propose several keystone ex- 
periments, including interference of superflows 
from two laser spots, vortex formation, scanning 
of moving dipole moments, and a giant increase 
of the luminescence distance by applying one- 
dimensional confinement potential. These exper- 
iments combined with our theory will open a new 
avenue for creating and controlling superflow of 
coherent excitons on nanoscale. 

The dramatic appearance of luminescence rings with 
radii of several hundred microns in quantum well struc- 
tures has been originally attributed to a boundary be- 
tween a positive hole gas diffusing from the laser spot 
and a negative electron gas located well outside the spot 
0, This implies an energetically unfavorable sepa- 
ration of the charges into hole-dominant positive and 
electron-rich negative regions on a macroscopic scale. 
The charge separation model reproduced some basic fea- 
tures of the luminescence ring formation. However, there 
are several observations that are not accounted for by this 
simple model including a non-monotonic dependence of 
the ring radius on laser power and significant short-range 
Coulomb interactions of the carriers [6]. The model is 
also refutable on the grounds of a rather short exciton 
formation time, typically ten picoseconds or less for rele- 
vant densities (> 10 10 cm -2 ) of photo-carriers @~E1- In- 



deed, photoelectrons (mass m) almost immediately emit 
optical phonons with the frequency lj, so that their speed 
is capped at about v = \j2fojj jm s» 4 x 10 7 cm/s. Thus, 
the exciton formation mean-free path turns out to be 
only a few microns, which is too short for creating the 
charge separation on the macroscopic scale of a few hun- 
dreds of micrometers far outside the spot. This is con- 
sistent with measurements of a photoluminescence (PL) 
spectrum in the magnetic field which have shown a dia- 
magnetic shift of the PL energy peak, indicating that the 
electrons and holes are bound into thermalized pairs (i.e 
excitons) in the region about a few tens of micrometers 
from the center of the spot 10]. As one gets closer to 
the center of the spot, this shift is changed from blue 
to a red one indicating an electron-hole plasma inside 
the spot, where electron- hole pairs rapidly form hot exci- 
tons by emission of phonons. Therefore, excitons, which 
can condense in a coherent superfluid, are formed from 
electron-hole plasma at rather short distances from the 
center of the laser spot, Fig.([T]). 

Importantly, since the electron-hole plasma expands 
inside the spot due to the Coulomb repulsion [ll| . en- 
ergy and momentum conservation requires that exci- 
tons initially form in states with finite center-of-mass 
momenta, K, which do not couple to light directly. 
If the lattice temperature, T is below the Berezinskii- 
Kosterlitz-Thouless (BKT) transition temperature, T< 
T^t w h 2 n ex /kBin ex , which is about a few Kelvin for 
the relevant exciton densities, n ex > 10 10 cm~ 2 , part of 
excitons should form a coherent state with a certain value 
of their momentum, K = Kq taken from electrons and 
holes forming pairs. Thus one can expect transitions of 
electrons and holes from their plasma into a flowing su- 
perfluid state of composed bosons in contrast to the more 
conventional stationary BEC of ultra cold atoms. As a 
result we expect an excitonic radial supercurrent appear- 
ing at some distance, i?o inside or just around the laser 
spot, J s (i?o) = Jo 7^ 0- This new state of matter allows 
excitons to travel on a macroscopic distance from the 
laser spot producing PL far away from the photo-excited 
region (Fig. 1), as observed in Refs.[l|, |2j. 



Model 

To describe the concomitant free superflow of two- 
dimensional (2D) Bosc-liquid of excitons we apply a 
generalized Gross-Pitaevskii equation [l3[ for the order- 
parameter, ip(R.) taking into account the exciton recom- 



FIG. 1: (Color online) Superfluid of slowly decaying excitons 
travels a macroscopic distance from the laser spot producing 
PL, when it stops far away from the photo-excited region. 
Relative length of vectors changes consistently with equation 
for the exciton momentum. The momentum k(p) drops 
to zero at the luminescence ring shown as a red ring. See also 
animation fl2|. where velocities in excitonic superfluid scale 
according to Eq.©; excitons start flashing in the animation 
as soon as their momentum drops below K c and they emit 
photons. 



bination rate 7/fi,, as originally proposed by Keldysh [l4| : 

- + V^(R) + F|V>(R)| 2 V(R) - = 0. 

(l) 

Here /x and V conveniently parameterize the average su- 
perfluid density and the short-range repulsion, respec- 
tively. 

Introducing dimensionless amplitude, /(r) = 
(F/^) 1 / 2 |V'(R)|, and the current density, j(r) = 
/(r) 2 V</>(R), Eq.([T]), is reduced to the following two 
equations 



and 



A/ ~ / + f + = 0, 



V-j = -/3/ 2 



(2) 



(3) 



where <fi(R) is the phase of the order parameter and coor- 
dinates r = R/£ are normalized on the coherence length 
£ = (h 2 /2m ex (i) 1 / 2 . Eq.@ is the familiar Ginzburg- 
Landau equation in the presence of the current and 
Eq. ([3]) is the continuity equation with the dimensionless 
decay rate, (i — 7/2/x. 

These equations are grossly simplified in the case of a 
low decay rate, f3 <C 1. In this case the amplitude of the 



FIG. 2: (Color online) Supercurrent density and exciton mo- 
mentum as a function of the order parameter. At relatively 
small exciton recombination rates the stable uniform super- 
flow (to the right from the vertical dotted line) is realized 
for lower momenta, k < l/\/3 and higher order parameters, 
f 2 > 2/3. The dot corresponds to the luminescence region, 
where excitons "burn" radiating light. Note that the order 
parameter (i.e. the superfluid density) increases approaching 
the luminescence region. 



order parameter and the current density change at the 
length scale of the order of 1//3 ^> 1 much larger than 
the coherence length, so that the first term ("kinetic pres- 
sure") in Eq.@ is negligible. Then the current density 
is expressed via the amplitude as 



(4) 



There are two brunches described with Eq.Q, but only 
one, which counter-intuitively corresponds to a larger or- 
der parameter f 2 2/3 and smaller momentum k = 
\J\ — f 2 l/\/3, is stable for small enough /3, Figj2](see 
supplementary material). 

A. One- dimensional (ID) superflow and PL. 

Let us first discuss the stable superflow in ID struc- 
tures which can be made by applying a ID confinement 
potential to experimentally available 2D samples. By us- 
ing the boundary condition, j(0) = V(m ex /2fi 3 ) 1 ^ 2 J , 
and assuming that Rq/£ <C 1//3, Eq. ([3]), written in the 
form dk/dx = /3(k 2 — 1)/(1 — 3k 2 ) for the exciton mo- 
mentum k(x), can be easily analytically solved: 



3(fc -k) + In 



(l + *c)(l-fco) 



(l-fe)(l + fco) 



= fa, 



(•5) 



where ko is found from k — k^ = j(0). 

To determine location of luminescence, we use that 
excitons can radiate in ID and 2D by resonant emission of 
photons only if their momentum is inside the photon cone 
[l5l [III , namely if superfluid center-of-mass momentum 
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is small enough: K ^ K c — E g y/e/c. Here E g is the 
semiconductor band gap and e is the dielectric constant. 
The exciton PL intensity is determined by the fraction 
of excitons inside the cone, where their decay rate /? is 
strongly enhanced [16| . For a narrow photon cone, K c <C 
Kq, the position of a bright luminescence stripe, x r i ng is 
found by taking k = in Eq.©, 



/r 1 



3fco — In 



1 + fep 
1-fco 



(0) 



v ring 



is macro- 



Since k is of the order of one, the distance x r 
scopic, i.e. large compared with the coherence length, as 
» 1. 

B. Controlling PL in 2D superflow. 

Let us now consider an already implemented experi- 
mental situation 0, Q , where excitons created around the 
central laser spot propagate freely in a plane. In contrast 
to the original proposal of Ref . [1| we suggest a continu- 
ous super&ow of excitons out of the excitation spot under 
steady-state photoexcitation, rather than their normal 
drift and diffusion. For isotropic 2D superflow the or- 
der parameter depends only on the distance p from the 
spot. Then 2D continuity equation (written for the 
exciton momentum k(p) = y/l — f 2 (p)) takes a simple 
one-dimensional form, 



dk 
dp 



k - k 3 



= 13- 



1 



p(l-3fc 2 ) ^l-3fc 2 



(7) 



The stable numerical solution of this equation is shown 
in Fig. [3] together with k(x) dependence of the momen- 
tum in ID structures using Eq.©. The analytical so- 
lution of linearized equation ([7]), k' p + k/p = —ft, which 
is k(p) — kopo/p— /3p/2 practically coincides with the 
numerical solution of the nonlinear Eq.([7]) in the whole 
relevant space (po is the distance from the center of the 
spot where the superflow starts). As in the ID case, 
PL ring radius corresponds to k(p r i ng ) = 0, so that 
Pring = (2fco/°o//?) ■ However, different from ID ge- 
ometry, where x r i ng scales as 1//3, the ring radius scales 
as pring oc 1 / \fj3 in the 2D case (see inset in Fig.((3|)). 
Therefore, the radius of PL for 2D case occurs to be 
much shorter than a distance between the laser spot and 
ID luminescence stripe. This prediction can be verified 
by channeling the superfluid, that is by applying ID con- 
finement potential in order to controllably transform 2D 
to ID geometry and check if PL shifts away from the 
spot. 

In dimensional units the PL ring radius is R r i ng ~ 
\J2YlRqKqI (m ex j). The distance from the center of the 
spot, Rq, where the superflow starts with the momentum 
Kq , can be estimated using the energy conservation. The 
average electron potential energy inside the spot is Vn, 
where V ~ 2irad, d is the effective separation between 
the electron and hole layers (« 10 nm), a = e 2 /(4ireoe) « 
2.3 x 10~ 28 /e Jm, and e is the static dielectric constant of 
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FIG. 3: (Color online) The exciton momentum as a function of 
the distance from the laser spot obtained using Eqs. H5I7|) for 
three different values of the decay rate /3 — 10 -2 , 10 -3 , 10 -4 
in both ID and 2D geometries. Dotted black lines correspond 
to solutions of the linearized 2D continuity equation. They 
practically coincide with the numerical solutions of the exact 
equation except of very short distances form the laser spot. 
Inset shows the dependence of the position of the lumines- 
cence region for these decay rates in log-log scale. It is seen 
that x r i n g scales as while pring scales as l/y^fl. This 



results in x r 



> Pr 



that is luminescence stripes formed 



in ID structures locate much farther from the excitation spot 
comparing with luminescence rings in 2D geometry. Gradu- 
ally changing the geometry from 2D to ID should move the 
luminescence farther from the laser spot. 



the semiconductor. It transforms into the exciton kinetic 
energy, i.e. Vn — K 2 /2m ex , so that Kq w V ' Anadnm ex . 
The radius Rq, where the superfluid is formed is about 
Rq = vr ex , where r ex is the exciton formation time. As 
a result we find 



R'T 



I hvT ex (16irdan / 'mex) 1 / 2 
7 



(8) 



Importantly 7 in this estimate is the nonradiative re- 
combination rate since the coherent excitons have their 
momentum outside the photon cone. The nonradiative 
lifetime, tq = ft/7 is an order of magnitude or more 
longer than the exciton radiative lifetime (l7| . Remark- 
ably, using the realistic material parameters T ex = lps 
@, tq = 10 ps [l7j], m ex = 0.2m e , e = 13 in Eq.© 
yields R r mg about 500 pm for the photoelectron den- 
sity n = 10 10 cm~ 2 , explaining the macroscopic radius 
of the rings. Dissipation processes stabilize the steady- 
state superflow of coherent excitons (see supplementary 
material), in contrast with a normal state flow, where 
individual excitons can scatter to a lower energy state 
emitting acoustic phonons. 

Also the dependence of the ring size and its PL intcn- 
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sity, I, on the excitation power, P, can be readily under- 
stood in the framework of the free superflow. According 

to Eq.ljHJ the ring radius scales as Rring oc n 1 / 4 -^/ 2 . The 
ring luminescence intensity, I is proportional to the ex- 
citon density on the ring, n ex . Excitons are pumped into 
the coherent state at i?o with the rate n/r ex , where the 
exciton formation time, T ex is inverse proportional to the 
photocarrier density n and it strongly depends on the 
exciton momentum K [9(. After flowing the distance 
Rring excitons decay radiating light, so their density on 
the ring scales with the photocarrier density 
Rqii/R 



■ring 1 ex • 



Parameterizing the momentum depen- 
dence of the exciton formation time as r ex oc Kq /n yields 
the following scaling: R r ing oc n^ 7 '^ 1 ^ 4 and I oc ■n/ b ~ r ^ 4 . 
Hence, we find R ring oc P^- 1 )/ 4 , J oc p( 5 - r )/ 4 in the 
case of a linear photoelectron population (n oc P), and 
Rr in g oc P^- 1 )/ 8 , I oc p( 5 - r )/ 8 in the case of n oc P 1 ' 2 . 
Since r can be a large number (^5) Q, our theory pre- 
dicts an increase of the ring radius and a decrease of the 
PL intensity with higher excitation power. In particular, 
the product R r i n gl is proportional to n for any expo- 
nent r independent of modeling the exciton-formation- 
time, thus resulting in R r ingl oc P 1 / 2 as observed [ij for 
n oc P 1 / 2 . More generally the dependence of r ex on the 
exciton momentum is essentially non- monotonous Q , so 
that the ring radius and the PL intensity may depend 
on the excitation power in a more complicated fashion as 
also observed @. 

There is a threshold density of photoelectrons n c in 
our theory in agreement with the experiments [J , and 
hence the laser power below which the ring does not 
form, which can be found from Kq = K c . At this and 
lower excitation power coherent exciton superfluid radi- 
ates just at the point of its formation inside or close to 
the laser spot. Using our estimate of Kq one obtains 
n c = P 2 e/(47rc 2 adm ex ) w 0.2 x 10 10 cm -2 . 

While a detailed description of the outer ring fragmen- 
tation into a periodic array observed in Ref. [Ij is outside 
the scope of the current paper, the periodic patterns are 
anticipated in our theory since the uniform coherent state 
of the outer ring is unstable due to a strongly enhanced 
recombination rate inside the photon cone (see supple- 
mentary material). 



C. Two-spot superflow and dark matter. 

Experiments with two rings created by spatially sepa- 
rated laser spots revealed that the rings attract one an- 
other, deform, and then open towards each other ||. 
This happens before the rings coalesce into a common 
oval-shaped ring, suggesting the existence of a "dark 
matter" outside the rings that mediates the interaction. 
Such a ring attraction is hard to consistently interpret in 
the classical electron- hole plasma model [TTJ] , where outer 
electron-rich regions from two spots naturally should re- 
pel each other. Here we show that this dark matter is 
the coherent superflow of excitons from two laser sports 
interfering with each other. Our continuity equation ([3]) 



becomes now an equation for the two components of the 
exciton momentum k = (k x , k v ), 



9(1 - k 2 )k x 3(1 - k 2 )k v 



dx 



dy 



= /?(fc 2 -l), 



(9) 



Using k = it can be reduced to a nonlinear partial 
differential equation for the phase of the order parameter 
4>{x, y) with the boundary condition k = fco on two small 
circles inside or just around every spot. As numerically 
shown above, the linearized continuity equation describes 
well the superflow from a single spot on the relevant dis- 
tance, Fig. ([3]), so that we consider a linearized version of 
Eq.©: 



A 2( j)(x, y) = -0, 



(10) 



where A 2 = d 2 /dx 2 + d 2 /dy 2 . The individual ring size 
is about ykopojp, so that one expects two independent 
PL rings when the distance d between the spots is large 
d ^> yj kopo/f3. The interaction should appear at the 
distance comparable to the ring size. If we assume that 
the superflow starts at a small distance from the centers 
of the spots such that po ^ &o/A then a solution of 



4>{x,y) = 



+ fc 0/ 9 ln(pip2) (11) 



approximately satisfies the boundary conditions for 
d/2 < yfk Q pl/(i. Here p h2 = y/(x ± d/2) 2 + y 2 is the 
distance from each spot, respectively. If we further as- 
sume that po < k 2 /((3ko), then the condition for PL, 
k ^ k c , is satisfied at the contour found from 



1 



1 



P 



{x-d/2) 2 +y 2 (x + d/2) 2 + y 2 2k oPo ' 



(12) 



and shown in FigfJJ The series of the obtained con- 
tours agree well with the experimental observations 
as the spots are brought closer, the rings deform, and 
then open towards each other. This happens well before 
the rings coalesce into a common oval-shaped ring (right 
lower panel in Fig.Q). 

Apart from the normal state diffusion in the charge- 
separation model, an alternative "superfluid" interpre- 
tation of PL rings in double and single QWs has been 
hampered by the lack of the theory of the free 2D super- 
flow. Based on our theory, we suggest that the macro- 
scopic luminescence rings and the dark matter in QWs 
originate in the slowly decaying superflow of excitons out 
of the laser spot at temperatures below their BKT tem- 
perature. The existence of a critical temperature about a 
few Kelvin of the macroscopic ring onset [5] strongly sup- 
ports our conclusion. Our results imply that Butov et al. 
[l| and Snoke et al. 0] discovered a new state of matter: 
the exciton superflow in the dark state. This state can 
be further experimentally verified with light scattering, 
STM or any other spectroscopy sensitive to the moving 
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dipole moments. Another piece of evidence for the super- 
flow should be an observation of a quantum interference 
and/or vortex structures with the characteristic coher- 
ence length £ = h(2m ex n ex V)^ 1 ^ 2 of about 15 nm. Last 
but not least prediction of the theory is a giant increase 
of the distance to the luminescence region by applying a 
ID confinement potential tuning the geometry from 2D 
to ID. 



FIG. 4: (Color online) PL contours from two laser spots 
for different distances d between the spots (measured in 
£\/kopo/P) obtained by using Eq. (|12[) . As in the experiment 
5] the rings first attract one another and, then, open towards 
each other. See also animations, where PL patterns from two 
laser spots separated by a fixed distance evolve with increas- 
ing k po//3 from 0.01 to 0.18 0. 
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Instabilities of the stationary superflow can be 
analyzed using a generalised time-dependent Gross- 
Pitaevskii equation, written in our units as 



-(A+l)^(r,t)+^ 2 (r,t)^*(r,i)-^V(r,t) = (*-u) 



dtp(r,t) 
dt ' 

Here u accounts for any kind of dissipation, for instance 
emission of plionons as in the time-dependent Ginzburg- 
Landau equation 

One can linearize this equation with respect to 
small fluctuations of the uniform superlow V's(r) = 
/(r) exp[i0(r] satisfying the stationary GP equation. 
Taking ip(x,t) = ip s (x) + ip(x,t) and ip(x,t) = 
r](x,t)i(} s (x) yields for small \rj(x,t) \ -C 1 

- V'L - IvUr'fL + i<P' x ) + f(v + V*) =(i~ u)f, (2) 

in ID geometry. The stationary amplitude f{x) and the 
phase 4>(x) satisfy the following pair of equations: 



-f" 

J XX 

f<t>'L 



f + f 
-pf, 







(3) 



where the second equation is the continuity equation. At 
small (J « 1 one can neglect the second derivative of 
the stationary amplitude f" x cx 1 , so that <j)' x = (1 — 
f 2 ) 1 ' 2 ee k and f x = -f3fk/(3k 2 -l), which substitution 
into Eq.© yields 



dk 



3fc 2 - 1 



+ (l-fc 2 ) (77 + 77*) = (i-u)f) 

(4) 

For a sufficiently long time, t < 1//3 ^> 1, a solution 
can be found in the form: rj(x,t) = aexp[z(cj/j — qx)] + 
&exp[— i(u>*t — qx)], where coefficients a and b satisfy the 
generalized Bogoliubov equations, 



(1 - k 2 - q 2 - 2qk 
(1 - k 2 -q 2 + 2qk 

with the eigenvalues, 



2if3qk 

3fc 2 - 1 
2if3qk 

3k 2 - 1 



+ uj + iuu)a = (k 2 — 1)6, 
uj + iuuj)b = (k 2 — l)a, 



2kq± W(l - k 



2i/3qk 
3fe 2 - 1 



iuu) 2 



(1-k 2 ) 2 . 
(5) 
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FIG. 1: Domains of stable and unstable uniform superflow. 
Increasing the exciton decay rate /3 relative to the dissipa- 
tion rate u shrinks the stability region to lower values of the 
coherent momentum k. 



In the absence of any decay and dissipation we recover 
the familiar Bogoliubov spectrum of excitations, shifted 
by the current, 



= 2kq ± vy + 2(l-fc 2 )g 2 



(6) 



The decay (3 and the dissipation u lead to an imaginary 
part of w, which is 



2/3kq \ (l-k 2 + q 2 ) 
3k 2 -l) ^qi + 2(1 - k 2 JP 



(7) 



for small u,j3 <C 1. The instability appears when 3w 
is negative. In the absence of the decay (j3 — 0) , this 
condition corresponds to the negative w^, or the positive 
LUg, which is the case if 



4fc 2 ^ q 2 + 2(1 - k 2 ). 



(8) 



Hence the uniform superflow with sufficiently small mo- 
mentum, k ^ 1/V3 is stable, while the higher momentum 
states with k > I/a/3 are unstable. A finite decay /3 ^ 
shrinks the region of the uniform superflow to k < 1/ 
as shown in Fig.l. 
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